Complex Impedance as a Diagnostic Tool for Characterizing Thermal Detectors 
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The complex ac impedance of a bolometer or microcalorimeter detector is easily measured and 
can be used to determine thermal time constants, thermal resistances, heat capacities, and sensitiv- 
ities. Accurately extracting this information requires an understanding of the electrical and thermal 
properties of both the detector and the measurement system. We show that this is a practical 
method for measuring parameters in detectors with moderately complex thermal systems. 
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I. INTRODUCTION 

Thermal detectors are used in a number of fields rang- 
ing from particle and plasma physics to astrophysics. The 
two commonly used forms of these detectors are bolome- 
ters to measure incident power and microcalorimeters to 
measure total energy. Despite these different applications 
their construction and operation principles are similar. 
In general, they are composed of an absorbing element 
to collect incident radiation or particles, a resistive ther- 
mometer coupled to the absorber, and a weak thermal 
link connecting the thermometer to a heat sink. 

The simple model treated in most analytical descrip- 
tions consists of a lumped heat capacity connected to a 
heat sink through a weak thermal link [Fig. 1(a)]. The 
theory predicting the responsivity, noise properties, and 
energy resolution of these detectors has been thoroughly 
developed. '^^^ However, more complex thermal models 
are often needed to describe real detectors. The inter- 
nal construction of these detectors often introduces ad- 
ditional time constants and noise terms which can have 
major detrimental effects on performance.^"* 

A detector's performance can be predicted and opti- 
mized for specific applications'*^*"^*' if the thermal and 
electrical properties of their components are known. 
Thermal conductivities, heat capacities, and sensitivi- 
ties can be extracted from measurements of the resis- 
tance vs. temperature relation, dc IV (current-voltage) 
curves, and time constants. ^''^'^^'^^ The thermal time con- 
stant can be found by measuring either the detector's 
impulse response function in the time domain^'^*'*'^'^'* 
or its complex ac impedance Z{uj) = dV/dl in the fre- 
quency domain. ^'^ In this paper we discuss the impedance 
method, which can often be used to find more precise pa- 
rameter values than the impulse response method.*^ 

We will show that the internal heat capacities and 
thermal conductivities of both simple and complex de- 
tectors, along with the thermometer temperature sensi- 
tivity and thermometer voltage (or current) sensitivity 
can be determined from measurements of dc IV-cuives 
and the ac impedance. As examples of this technique 
we fit impedance data to the ideal detector model and 
one of the more complex thermal circuits presented by 
Galeazzi and McCammon^. This complex model in- 
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FIG. 1: a) The ideal model of a thermal detector consists of a 
lumped heat capacity C connectecd to a heat sink through a 
weak thermal link G. b) A more realistic model that includes 
internal couplings." 



eludes electron-phonon decoupling*^"^* in the thermome- 
ter and a thermal resistance between the absorber and 
thermometer. 

We begin with a review of the impedance of the sim- 
plest detector model in Sec. II. Section III presents 
example fits of simple and complex thermal models to 
impedance data. In Sec. IV we discuss the effects of 
stray electrical capacitances and inductances in the bias 
and readout circuits. While the examples presented here 
utilize a voltage readout circuit, the equations can readily 
be transformed to current readout. 



II. THE SIMPLE THERMAL DETECTOR 

The simplest, or ideal, detector consists of a single 
lumped heat capacity C connected to a heat sink through 
a weak thermal link of conductivity G [Fig. 1(a)]. This 
model has been widely discussed by other authors;^""^'^ 
we briefly review their results here. 

The dynamic impedance for this simple model is given 

by 
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detector resistance 
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Joule power dissipated in detector = V"^ / R (3) 
thermal conductivity to heat bath = dP/ dT (4) 
thermal time constant = C/G (5) 
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For a linear thermometer the resistance is dependent only 
on its temperature and (iy = 0. However, this is not 
generally the case so we will retain the term here. 
Equation (1) can be rewritten as 
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Equation (8) describes a semi-circle in the complex plane 
of radius ||Zoo — Zq] centered on the real axis at ^{Zoc, + 
Zq). The frequency at the peak of the circle is given by 
the dynamic time constant fz- This time constant is not 
the same as the effective time constant, which describes 
the effect of electro-thermal feedback on the detector re- 
sponse and is dependent on the relative resistance val- 
ues of the detector and load resistor. The dynamic time 
constant is a property of the detector only. Equations 
(8)-(ll) are equivalent to those given by Mather^. 

When > [as is the case for a superconduct- 
ing transition-edge sensor (TES)] it is possible for the 
impedance to become infinite [Z/„ = 1 — /3„ in Eq. (10)]. 
To avoid this, one might instead measure the complex 
admittance A{u}) = 1/Z{u>). The relations describing the 
admittance are easily found by transforming the preced- 
ing relations using the dual circuit theorem, namely 
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C (parallel capacitance) 



V, (12) 
S{=l/R), (13) 
A. (14) 
L (series inductancejjl5) 



Since it is customary to keep a and (3 as derivatives of 
resistance rather than conductance, the signs of these 
quantities will change: 
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The complex admittance is then 
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With these equations we see that the path a negative 
detector {a^ < 0) traces through the complex impedance 
plane is equivalent to the path traced through the com- 
plex admittance plane for a positive detector (a^ > 0). 



III. IMPEDANCE MEASUREMENTS 

Figure 2 shows examples of impedance measurements 
and fitted models for three different microcalorimeter de- 
tectors. The first example is for a doped silicon thermis- 
tor with no absorber [Figs. 2(a) and 2(b)]. The single 
semi-circle is adequately represented by Eq. (8) with the 
parameters shown in Table I. The mismatch between 
the data and models at high frequency is most likely due 
to incomplete modeling of stray capacitances (see Sec. 
IV). In the other two examples, composite absorbers 
have been glued to the detectors [Figs. 2(c)-2(f)]. The 
more complex behavior of Z{lli) indicates significant ther- 
mal resistances which can be reasonably fit by the ther- 
mal model of Fig. 1(b). The complex impedance for this 
model is derived by Galeazzi and McCammon^ as their 
Eq. (110) (note that they use f3i rather than /3„). 

The complete characterization of these thermal detec- 
tors includes the sensitivities of the thermistor a and /3, 
and the heat capacities and thermal conductivities of all 
constituent parts. However, these parameters cannot be 
determined from measurements of the impedance alone. 
The impedance of the ideal thermal model is completely 
described by the three parameters Zq, Z^c, and in Eq. 
(8), which are determined by fitting the measured Z{uj). 
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FIG. 2: Examples of impedance measurements for a bare thermometer [(a) and (b)] and thermometers with absorbers [(c) (f)]. 
Data are shown as dotted lines and the fitted model as solid lines. The figures on the left [(a), (c), and (e)] show the frequency 
variation of the real and imaginary parts of the detectors' complex impedance while the figures on the right [(b), (d), and (f)] 
plot the same data in the complex plane. The divergence of the fits at high frequency is due to limitations of modeling the 
electrical circuit stray capacitances. The x:y axis ratio is 1:1 in (b), and 2:1 in (d) and (f). For clarity, data above 5kHz have 
been removed in (a) and (b). The fit model parameters for (a)-(d) are given in Table I. 
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TABLE I: Example detector parameters 



Parameter 


Bare" 


Absorber' 


Method" 


R [MO.) 


36 


1.3 


dc measurement 


P (pW) 


0.22 


7.0 


dc measurement 


T (mK) 


70 


140 


R-vs-T calibration 




-7.5 


-4.5 


R-vs-T calibration 


G (pW/K) 


35 


230"^ 


dc IV-curve 


Cp (pJ/K) 







fixed 


fz (ms) 


1.5 




Impedance fit 




3.3 


0.053 


Impedance fit 


Zoo (Mfi) 


34 


1.1 


Impedance fit 




-0.058 


-0.17 


Eq. (9) 


Lv 


-0.83 


-0.91 


Eq. (10) 


T (ms) 


2.7 




Eq. (11) 




-9.0 


-4.2 


Eq. (2) 


Ce (pJ/K) 


0.093 


0.13 


Eq. (5), Impedance fit*^ 


Ca (pJ/K) 




0.64 


Impedance fit 


(pW/K) 




460 


Impedance fit 



"Bare thermistor shown in Figs. 2(a)-(b) 
''Thermistor with absorber shown in Figs. 2(c)— (d) 
'^Method used to measure the indicated parameter 
''Series combination of G and Gg 

"^Ce for the bare thermistor is found from Eq. (5); it is fit directly 
for the thermistor with absorber. 



A dc determination of the resistance yields L^, and 
T through Eqs. (9) (11). However, C, and G cannot 
be separated from Ly and r without another independent 
measurement. For this last measurement, G{T) is fitted 

to dc IV -curves. 

The data in Figs. 2(c)-2(f) were fit using the model 
in Fig. 1(b), which separates the detector heat capacity 
into contributions from the absorber Ca, thermistor elec- 
trons Ce, and thermistor phonons Cp, connected by the 
labeled thermal conductivities. A complete characteri- 
zation requires eight parameters: three heat capacities, 
three thermal conductivities, and two thermometer sen- 
sitivities. In principle, this model contains three separate 
thermal time constants which could appear as three sep- 
arate circles in the complex impedance data. In practice, 
the phonon heat capacity in our doped silicon thermistors 
is much smaller than either the electron or absorber heat 
capacities, placing this third time constant at frequen- 
cies beyond the practical measurement range. Therefore, 
Cp « and G and can be replaced by their series com- 
bination. The remaining six parameters are still perfectly 
correlated in fits to the impedance alone so independent 
measurements are required to separate them, just as in 
the ideal model. 

The procedures oiitlined above for characterizing ther- 
mal detectors contain at least one known systematic er- 
ror. The fits of G{T) to /^-curves assume that = 0, 
while the impedance fits and resistance measurements 
clearly indicate that Py (see Table I). This prob- 



lem could be alleviated using an iterative procedure (re- 
fitting the /V-curves with the non-zero (3y) or by per- 
forming simultaneous fits to both the IV-ciivyca and 
impedance. Without these corrections the small val- 
ues of py measured in the doped silicon thermometers 
{\f3y \ < 0.2) introduce only a few percent uncertainty in 
the IV-cuive determination of G{T). These required cor- 
rections may be larger in TES detectors which can have 
relatively large values of Pi (~ 1-3).^^ 



IV. READOUT CIRCUITS 

A. The Transfer Function 

Figure 3 illustrates two equivalent circuits for biasing 
thermal detectors, either of which can be used to measure 
the impedance. The first utilizes voltage readout [Fig. 
3(a)], the other current readout [Fig. 3(b)]. Our tech- 
nique for measuring the dynamic impedance is to add a 
small ac signal to the dc bias and measure the resulting 
complex transfer function as a function of frequency. The 
transfer functions are defined as T(w) = Kut(w)/Vin(w) 
for voltage readout and T{lo) = /out('^)/-^in(w) for cur- 
rent readout. The ac signal can be either a random noise 
source with a bandwidth spanning the frequencies of in- 
terest or a sinusoidal source that can be scanned through 
the desired frequency range. We use a commercial spec- 
trum analyzer to measure the complex transfer function 
using both of these methods. The sine-sweep method 
could also be implemented with a two-phase lock-in am- 
plifier and swept sine generator. 



B. Electrical Strays 

For an ideal bias circuit the transfer function would 
simply be the voltage or current divider formed hj Rl 

and Z{lu). In practice, additional reactive components 
exist in the form of stray capacitances and/or induc- 
tances. These stray reactances will introduce additional 
phase and amplitude shifts in T{lo), but it is possible to 
correct for these effects if the strays in the bias circuit 
can be accurately modeled and measured. 

The stray reactances can be determined by measuring 
the transfer function of the circuit with the reactive part 
of Z(ijj) removed. This is done by either replacing the 
detector with a pure resistor or by measuring the trans- 
fer function of an unbiased ((Vin) = (/in) = 0) detector. 
If the unbiasc;d transfer function is measured with a suf- 
ficiently small ac signal, such that the detector dissipates 
negligible power (P » 0), then the detector is purely 
resistive [see Eqs. (1) and (2)]. Any observed reactive 
component under these conditions must be due to stray 
reactances in the circuit. 

The dominant reactances in our voltage readout cir- 
cuit for silicon detectors are shunt capacitances across 
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FIG. 3: Bias circuits for measuring detector transfer functions. Both circuits consist of a thermometer impedance Z{ui), 
resistance R, and load resistor Rl- a) Voltage readout circuit for doped silicon thermometers. The capacitors connected by 

dotted lines are used to model stray capacitances in the system. Cl and Ct are shunt capacitances across the load resistor and 
thermometer, respectively. The load resistor is one physical resistor, but it has been split in this schematic {Rl = Rli + RL2) 
in order to partially model the distributed capacitances Ch- b) Current readout circuit for TES. The inductors shown as dotted 
lines axe used to model stray inductances in series with the thermometer, Lt, and load resistor, Ll. If Cb = these two circuits 
are dual equivalents and all equations describing one system can be transformed to the other system using the substitutions in 
Eqs. (12)-(15). 



the detector and load resistor (Ct and Cl in Fig. 3, re- 
spectively). For this model circuit the transfer function 
is given by 



T(c.) 



Z{oj) l+jWTL 



Mn Rl + Z{uj) l+jWTt 



where 



Th = 

n = 



RlCl and 



(23) 



(24) 
(25) 



This simple model for the electrical strays also traces 
a semi-circle in the complex plane if the detector 
impedance is completely real, Z{u!) — R. 



T(a;) = ^°±^ + ^ 



T^. 1 



- 3^n 



where 



lim T{uj) = 



-JWTt 



and 



To = T(0) = 



R 



-Rl + R 



(26) 

(27) 

(28) 



The dual circuit theorem can be used to transform 

these equations from the voltage readout circuit of Fig. 
3(a) to the current readout circuit of Fig. 3(b). Equations 



(23)-(25) are then 
T{oj) = 

Jin 



A{u;) 



1 



-JWTL 



Si^ + A{uj) l-FjwTt 
-Rl 1 -I- jg^TL 



where 



TL = -Z^l/-Rl and 



n = 



Rl + Z{uj) ■ 



(29) 

(30) 
(31) 



Equations (26)-(28) can be similarly transformed. 

Real strays are complex, involving distributed reac- 
tances rather than (or in conjimction with) the simple 
parallel capacitances and series inductances discussed 
here. It is the accuracy of the circuit model, not the 
values of the circuit elements (both real and stray) that 
limits the maximum useful frequency for impedance mea- 
surements. This will particularly limit the accuracy with 
which Zoo and Py (or Aoo and /?») can be measured. 

No stray inductances have been included in Fig. 3(a) 
as they are generally negligible for high-impedance sili- 
con detectors. Similarly, the effect of stray capacitances 
can be neglected for the very low impedance TES de- 
tectors. In this case the stray inductances shown in Fig. 
3(b) can be significant. However, these stray inductances 
still seem to be less of a problem for the low-impedance 
detectors than stray capacitances are for high-impedance 
detectors. As a result, the impedance of a TES can often 
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FIG. 4: Example of an unbiased transfer function, a) The 
measured complex transfer function (diamonds represent ev- 
ery tenth data point) and fits to two stray capacitance models 
[dotted and dashed lines; Fig. 3(a)]. If no strays were present 
the data would follow the solid lines, b) The inferred detec- 
tor impedance after correcting for each stray model (dotted 
and dashed lines) is compared with an ideal resistor (solid 
lines). One model uses only shunt capacitances across the de- 
tector and load resistor (Cb ~ 0) while the other also includes 
a distributed capacitance within the load resistor (Cb 7^ 0). 
The impedance data above 5 kHz have been removed from the 
bottom plot for clarity. 

be measured up to higher frequencies (tens of kHz)^^'^^ 
than the silicon detectors (a few kHz). 

Wc note that the utility of the impedance method is 
not limited to only those systems which can be described 
by the bias circuit models in Figure 3. It is easily ex- 
tended to any bias circuit whose frequency dependent 
transfer function can be modeled and measured. This 
includes bridge circuits and long transmission lines, such 
as those often employed in fusion research. ^'^^'^^'^^ 



C. Measuring Stray Capacitances 

The shunt capacitance across the silicon detectors is 

typically Ct ^ 10 pF. The load resistors used in this work 
are made of nickel-chromium thin film deposited on a 



FIG. 5: Measured and fitted transfer functions of a biased 

detector with absorber, a) Transfer function, b) Detector 
impedance inferred from the transfer function using the mea- 
sured stray capacitances. 



thin silicon nitride insulating layer on a degenerate silicon 

substrate^'^. These were found to have shunt capacitances 
on the order of Cl ~ 1 pF. A set of Nichrome resistors 
on quartz substrates^'^ had shunt capacitances an order 
of magnitude smaller. 

These two shunt capacitances alone are not always suf- 
ficient to fit the unbiased transfer functions at high fre- 
quencies. In an effort to better model the distributed ca- 
pacitance in the load resistor the resistance is split into 
two pieces and a third shunt capacitor is added, Cb [Fig. 
3(a)]. Typical values of these additional components are 
Ch < IpF and RlJRl2 « 1 - 5. 

Example fits of an unbiased transfer function are 
shown in Fig. 4. Including the effects of the stray ca- 
pacitances in fitting the unbiased data should result in a 
completely real (resistive) impedance which is indepen- 
dent of frequency. The lower panel of Fig. 4 makes it clear 
that our simplest stray model (Cb = 0) accomplishes this 
only at frequencies below 500 Hz. If we include the dis- 
tributed capacitance in the load resistor (Cb ^ 0) the 
effects can be removed up to frequencies of » 1 kHz. The 
remaining deviations above 1 kHz indicate that our stray 
models are not an adequate representation of the real 
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circuit at these frequencies. 

D. Biased Transfer Functions 

In Fig. 5(a) we measure the transfer function of a bi- 
ased detector (dotted hnes). Using this transfer function 
and the measured values of the load resistor and stray ca- 
pacitances wc can calculate the impedance of the thermal 
detector from Eq. (23) or its equivalent; this impedance 
is shown as the dotted lines in Fig. 5(b). The best fit 
to the biased transfer function [solid line in Fig. 5(a)] is 
obtained by varying the detector parameters (i.e. a, /3, 
C"s, G's) while the strays remain fixed at their measured 
values. The impedance for these best fit detector param- 
eters is shown as the solid lines in Fig. 5(b). As already 
mentioned in Sec. Ill, the divergence at high frequencies 



is most likely due to inadequate modeling of the stray 
capacitances. 
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